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1. Introduction 



The aim of lattice QCD calculation is to nonperturbatively evaluate physical quantities from 
first principles with the systematic errors under control. Most troublesome are quenching effects 
and uncertainties associated with chiral extrapolation, which should be removed by 2+1 flavor 
simulation at the physical point. 

The previous CP-PACS/JLQCD project^] focuses on performing Nf = 2+1 lattice QCD 
simulations incorporating the dynamical strange quark by the Polynomial Hybrid Monte Carlo 
(PHMC) algorithm^], where the nonperturbatively 0(a) -improved Wilson quark action[||] and the 
Iwasaki gauge action[|[] are employed on a (2 fm) 3 lattice with the up-down quark mass down to 
67 MeV. This project is succeeded by the PACS-CS (Parallel Array Computer System for Com- 
putational Sciences) project^, ^, [7[ [| [| [l0|, [jl], [l2|, [B, 14j, which aims at the physical point 
simulation with the same quark and gauge actions as the CP-PACS/JLQCD project on enlarged 
physical volumes. We reduce the up-down quark masses using the domain-decomposed HMC 
(DDHMC) algorithm armored with several improvements and implement the UV-filtered PHMC 
algorithm^] for the strange quark. 

In this report we demonstrate the feasibility of a direct simulation at the physical point and 
examine viability of extrapolation method with the chiral perturbation theories. We apply the 
SU(2) and SU(3) chiral perturbation theories (ChPTs) to the pseudoscalar meson masses and the 
decay constants with the up-down quark mass ranging from 3 MeV to 24 MeV and two choices 
of the strange quark mass around the physical value. The low energy constants in the SU(2) and 
SU(3) ChPTs are determined and compared with those obtained with the different quark and gauge 
actions employed by other groups. We discuss the convergence behaviors of the SU(2) and SU(3) 
chiral expansions up to NLO. We also investigate the quark mass dependence of the nucleon mass 
employing the SU(2) heavy baryon chiral perturbation theory up to NNLO. We finally make a 
comparison between the physical quantities directly measured on the physical point and those 
estimated by the extrapolation method. This comparison reveals the necessity and the effectiveness 
of the direct simulation at the physical point. 

All the calculations are done using the PACS-CS computer with a total peak speed of 14.3 
TFLOPS developed and installed at University of Tsukuaba on 1 July 2006. A part of the results 
are already reported in Ref. [12, 13, 14]. 



2. Simulation details 



Our simulation is carried out using the 0(a) -improved Wilson quark action with a nonpertur- 
bative improvement coefficient csw = 1 .715 [^J] on a 32 3 x 64 lattice. The lattice spacing is found 
to be 0.0907(13) fm whose determination is explained later. Simulation parameters are summa- 
rized in Table |]. We choose seven combinations of the hopping parameters (k" U( j, K" s ) based on the 
analysis of the previous CP-PACS/JLQCD results. We employ two choices of K s to investigate 
the strange quark mass dependences. The heaviest combination (fc u d,fc s ) = (0.13700,0.13640) 
corresponds to the lightest one in the previous CP-PACS/JLQCD project. We expect that the com- 
bination (K"ud)Ks) = (0.137785,0.13660) would be the physical point, which is estimated at the 
early stage of our analyses. 
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Table 1: Simulation parameters. MD time is the number of trajectories multiplied by the trajectory length 
X. 
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The DDHMC algorithm! 15] allows us to access the small up-down quark mass region closer to 
the physical point. It is employed for the simulation points with fc u d < 0.13770. We divide the full 
lattice into 8 4 blocks which are used for a geometric separation of the up-down quark determinant 
into the UV and the IR parts. This separation makes possible to incorporate the multiple time scale 
integration scheme[|fj]], which reduce the stimulation cost significantly. We employ the replay 
trick[[TJ, 17 ] choosing the threshold AH > 2. The detailed description about the DDHMC algorithm 
is given m Refs. [fH, 0]. 

Reducing the up-down quark mass, however, the IR force Fir becomes less stable yielding 
spike-like fluctuations, which results in larger AH with higher replay rates. At the simulation 
point of (fC u d,K" s ) = (0.13781,013640) we incorporate the mass preconditioning! 18, [OJ to tame 
the fluctuations of the IR force, which is divided into the preconditioner F/ R with a new hopping 
parameter K^ d and the preconditioned part Fir. K^ d is parametrized as K^ d = pifc u d with Pi l ess 
than unity. The step sizes for Fq, Fjv> F/ r , Fir are controlled by four integers (No,Ni,N2,N^) as 
<5t g = t/(N NiN2Ns), <5ttjv = t/(NiN 2 N 3 ), 8t{ r = t/^A^), 8t{ r = t/Nt, with % the trajectory 
length. This algorithm is refereed to as MPDDHMC. We choose 5t s = 5t ir for the strange quark 
force in the UVPHMC algorithm based on our observation that | |F S | | w | |F IR | |. 

At (fC u d,K" s ) = (0.13781,013640) we also implement several improvements for the inversion 
of the Wilson-Dirac operator on the full lattice. First one is the chronological guess for the ini- 
tial solutions^] with the use of the last 16 solutions. Second one is a nested BiCGStab solver 
consisting of the outer solver with double precision arithmetic and the inner one operated with 
single precision. The latter with an automatic stopping condition from 10~ 3 to 10~ 6 works as a 
preconditioner for the former. We employ a stringent tolerance \Dx — b\/\b\ < 10~ 14 for the outer 
solver to retain the reversibility of the molecular dynamics trajectories to high precision. Third one 
is the GCRO-DR (Generalized Conjugate Residual with implicit inner Orthogonalization and De- 
flated Restarting) algorithm^]]] which is robust against the small eigen values of the Wilson-Dirac 
operator. It take over the inversion once the nested BiCGStab solver becomes stagnant. 

For the run at (fc ud , fc s ) = (0.137785,013660) we apply twofold mass preconditioning (MP2DDHMC) 
to the IR force Fir, which is decomposed into F IR , F IR and Ftr. Two hopping parameters are addi- 
tionally introduced: K^ d = P2K^ ld = p2PiK u d with pi and p 2 less than unity. We need five integers 
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(No,N\,N2,Ni,N4.) to adjust the step sizes for Fq, F\j\, F ir , F{ r , Fir. Our choice is <5t s = <5t ir for 
the UVPHMC algorithm. 

3. ChPT analyses on the pseudoscalar meson sector 

We examine the chiral behaviors of the pseudoscalar meson masses and the decay constants 
based on the SU(3) and SU(2) ChPTs up to NLO. Since the redefinition of some of the low energy 
constants (LECs) makes the one-loop expressions of the Wilson ChPT in terms of the AWI quark 



masses[|22|] equivalent to those in the continuumQ12|], we focus on the analyses with the use of the 
continuum ChPTs. 

3.1 SU(3) ChPT 

The one-loop expressions in the continuum SU(3) ChPT are given by[[23|] 

2 



^-=B {l+^-^ + ^ ( 16m ud (2L 8 -Ls) + 16(2m ud + m s ) (2L 6 - L 4 )) J (3. 1) 

Bo\ l + -^ + —^(S(m ud + m s )(2L s -L 5 ) + l6(2m ad + m s )(2L 6 -L 4 )) 1(3.2) 
(m ud + m s ) { 3 /o J 



2m ud 
nr. 



fx = fo\l- -Hk + ^ (8m ud L 5 + 8(2m ud + m s )L 4 ) \ , (3.3) 



[3 3 3 2B() 1 
fx = fo 1 1 - - - 4 Mj? + -j2 ( 4 ( mud + ms ) Ls + 8 ( 2m ud + I , (3-4) 



where we have six unknown LECs Bo,/o> ^4,5,6,8- Mps denotes the chiral logarithm defined by 

16^ \ n 

where 



2 

ml = 2m nd B , m\ = (m ud + m s )5 , m n = - (m ud + 2m s )B (3.6) 
with jU the renormalization scale. The pseudoscalar meson decay constants are calculated with 



the nonperturbatively 0(a) -improved axial vector current[24], though the renormalization factor is 



perturbatively evaluated up to one-loop level [25 ,26]. We determine the LECs by a simultaneous fit 



of m|/(2m ud ), m\/(m ud +m s ), f % and fx including the finite size corrections at one-loop level[^7|]. 
In Table g we compare our results for the LECs with the phenomenological estimates with 



experimental inputs [|27|,|28|] and recent 2+1 flavor lattice QCD results [|29|,|30|]. We observe that the 
situation is rather complex: Some LECs are consistent and others are not. To make the comparison 
easier we convert the SU(3) LECs to the SU(2) ones, where the number of LECs are reduced from 
six to four: B, f, l 3 , l A . We obtain J 3 = 3.47(11), l A = 4.21(11) and f 3 = 3.50(11), k = 4.22(10) 
with and without the finite size corrections, respectively. These results are plotted in Fig. [T] together 
with the phenomenological estimates and the recent 2 and 2+1 flavor lattice QCD results. We find 
that all the results reside in 3.0</3<3.5 and 4.0</4<4.5 except the MILC result for f 3 which shows 
exceptionally small value. 
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Table 2: Results for the LECs in the SU(3) ChPT together with the phenomenological estimates Q27|, [280 and 
the RBC/UKQCD p^] and the MILC results [|3C|]. /o is perturbatively renormalized at one-loop level. L4.5.6.8 
are in units of 10~ 3 at the scale of 770MeV. 



PACS-CS phenomenology |g7|, g§ RBC/UKQCD []29j MILC[Q 

w/oFSE w/FSE 



/otGeV] 


0.1160(88) 


0.1185(90) 


0.115 


hlh 


1.159(57) 


1.145(56) 


1.139 


U 


-0.04(10) 


-0.06(10) 


0.00(80) 


U 


1.43(7) 


1.45(7) 


1.46(10) 


1L(, — L4 


0.10(2) 


0.10(2) 


0.0(1.0) 


2L 8 -L 5 


-0.21(3) 


-0.21(3) 


0.54(43) 


% 2 /dof 


4.2(2.7) 


4.4(2.8) 





0.0935(73) 

1.33(7) 
0.139(80) 
0.872(99) 
-0.001(42) 
0.243(45) 
07 



1.21(5) (±3*) 
0.1(3) {tl) 
1.4(2) (+?) 
0.3(1) (jj) 
0.3(1)(1) 



1 




1 

KD) 




1 1 1 1 1 

PACS-CS(SU2) 






1 


HCH 


1 


1 1 

PACS-CS(SU2) 










PACS-CS(SU2-FSE) 








HCH 




PACS-CS(SU2-FSE) 






• 




PACS-CS(SIB) 








m 




PACS-CS(SU3) 






• 




PACS-CS(SU3-FSE) 








m 




PACS-CS(SID-FSE) 






HCH 
H»H 




RBC/UKQCD(SU2) 
RBC/UKQCD(SU3) 








• 




RBC/UKQCD(SU2) 
RBC/UKQCD(SU3) 










MILC(SU3) 








1— • — 1 




MILC(SU3) 


1 


• 


1 






















JLQCD 












JLQCD 






1 ■ — 1 










1 ■ — 1 








\-m-t 




ETM 








H 




ETM 






hhh 




CERN 








HH 




CGL 










GL 














GL 

1 1 


1 




• 

1 


1 






h 

1 


• 1 


1 


1 


1 


2 3 4 


5 


6 7 8 


2 


3 


4 5 


6 


7 



l 3 -bar 



1 -bar 

4 



Figure 1: Comparison of the results for h and I4. Black symbols denote the phenomenological estimates [ftlj 
|32] ]. Blue ones are for 2 flavor results [f33[ |34] , [35| l. Red closed (open) symbols represent the results for the 
SU(3) (SU(2)) ChPT analyses in 2+1 flavor dynamical simulations [0191. 



Although we obtain reasonable values for I3 and I4 in the SU(3) ChPT fit, there exists a 
concern that the values of # 2 /dof are unacceptably large. Its origin is found in the fit of f n and 
fx depicted in Fig. ^J, where the strange quark mass dependence between the data at (fc u d, fc s ) = 
(0.13754,0.13640) and (0.13754,0.13660) is not properly described by the SU(3) ChPT. The con- 
vergence of the SU(3) ChPT is checked in Fig. |3[ For fa the NLO contribution is about 40% of the 
LO one almost independent of the up-down quark mass. This implies that the strange quark mass 
is not small enough to be treated in the SU(3) ChPT up to NLO. However, the extension to NNLO 
is not a good choice in a practical sense: The NNLO contribution could be 20% by naive power 
counting and a full determination of additional LECs is not achieved without increasing the data 
points significantly. We take a different direction: the SU(2) ChPT fit with the analytical expansion 
in terms of the strange quark mass about the physical value. 
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Figure 2: Fit results for f n (left) and /a- (right). Black symbols represent the lattice results. Red triangles 
denote the SU(3) ChPT results plotted at the measured quark masses. Open and filled symbols distinguish 
between the results at ?c s = 0.13660 and 0.13640. Star symbols represent the extrapolated values at the 
physical point which is denoted by vertical dotted line. 
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Figure 3: Ratio of the NLO contribution to the LO one in the chiral expansions of the pseudoscalar meson 
decay constants. Solid (open) lines represent the SU(3) (SU(2)) ChPT case. The strange quark mass is fixed 
at the physical value. 



3.2 SU(2) ChPT 

We try further chiral analyses on the pseudoscalar meson sector employing the SU(2) ChPT 
up to NLO, where the K meson is treated as a matter field in the isospin 1/2 linear representation 



of the SU(2) chiral transformation[36]. The strange quark contributions to the SU(2) LECs are 
analytically expanded about the physical value. The SU(2) ChPT formulae for m % and f n are given 
by 
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with m n = m ud B. B and / are linearly expanded in terms of the strange quark mass: 
m^p and / = fs + m s fs . For niK and fa we employ the following fit formulae: 



m K = m K + Annuel, fK = f i 1 + Pf"hxd ~ J 
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Figure 4: Same as Fig. 2 for the SU(2) ChPT. 



A simultaneous fit to m 2 r /2m ut i, f n and fg leads to the results of It, = 3.14(23), I4 = 4.04(19) 
and Z3 = 3.23(21), I4 = 4.10(20) with and without the finite size corrections, respectively. We 
apply an independent fit to m\. The results for Z3 and l\ are added in Fig. |l|, where we find that 
all the lattice results are consistent except the MILC result of /3. The fit results for f n and 
are depicted in Fig. Q The lattice data are well described by the SU(2) ChPT combined with the 
analytic expansion in terms of the strange quark mass. The corresponding # 2 /dof are 0.43(77) and 
0.33(68) with and without the finite size corrections, respectively, which are reasonable magnitude 
contrary to the SU(3) case. In Fig. |3| we compare the convergence behaviors between the SU(3) and 
the SU(2) cases. The ratio of the NLO contribution to the LO one in the SU(2) ChPT is drastically 
reduced from the SU(3) case in the range of ra ut i<0.01. 



4. ChPT analyses on the baryon sector 

It is intriguing to investigate the quark mass dependence of the baryon mass based on the 
heavy baryon ChPT (HBChPT). Let us consider the nucleon mass as an example. The SU(2) heavy 
baryon ChPT formula up to NNLO is given by[^7|] 

A 2 6 8a 3 
mjv = mo — 4cim K — %~Fy m % 
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with mo,f Jl ,ci,C2,C3,gA,ei low energy constants. Note that the inclusion of the 0{m 5 n ) term, which 
is obtained by the relativistic baryon ChPT[38], little affects on the following analyses because of 
its tiny contribution. 

The HBChPT fit is performed under several constraints. Firstly, the fit parameters are re- 
stricted to mo, ci, e\. For a direct comparison with the 2 flavor results of the QCDSF-UKQCD 
Collaborations!!!], we choosey = 1.267,c 2 = 3.2GeV _1 wifhc 3 = -3.4GeV _1 orc 3 = -4.7GeV" 
The former is refereed to as Fit- A and the latter as Fit-B. We employ the value of f K at the chiral 
limit obtained by the SU(2) ChPT fit for m n ,f n and fa- Secondly, the ChPT fit is applied to the 
data series at fc s = 0.13640 which show a few % level of variation for m^ WI . We take two fit ranges: 
0.13781 < fc ud < 0.13727 for Range-I and 0.13781 < fc ud < 0.13700 for Range-II. 

Table [| presents the fit results for the LECs in comparison with the QCDSF/UKQCD results^] 
We also list the ETM results obtained at /3 = 3.9, which yields similar lattice spacing to ours, 
with the choice of c 3 = — 3.45GeV~'[^(|]. The results show reasonable agreement in case that 
we take the similar values for c 3 , We also give the nucleon sigma term evaluated through o^x = 
m 2 K (dniN / dm 2 ^) with the aid of Eq.(4.1). In Fig. || we draw the fit results and investigate each con- 
tribution of the LO, NLO, NNLO and terms in Eq.(|Q|). It is remarkable that the lattice 
results are well described up to m\ ~ 0.5 GeV 2 both for Fit- A and Fit-B. We also find reasonable 
values for # 2 /dof in Table [| However, this is achieved by a drastic cancellation between the LO 
and NLO contributions. It is hard to say that the convergence of the SU(2) HBChPT is controlled 
beyond m 2 ~ 0.2 GeV 2 . 




Figure 5: Fit results of the nucleon mass using the SU(2) heavy baryon ChPT formula (left) and each 
contribution of the LO, NLO, NNLO and 0(m 5 K ) terms in the case of Fit- A with Range-I (right). 
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Table 3: Fit results for niQ,c\,e\ in comparison with the QCDSF/UKQCD and the ETM results. Values for 
the nucleon sigma term are also presented. 





Range-I 


Rarrj 
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Fit-A 
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-3.4 
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-3.4 
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-1.00(10) 
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2.9(4) 


2.4(3) 
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0.1(0.9) 


0.0(0.5) 


0.3(0.9) 


1.2(1.6) 


o N7l [MeV] 


51.4(7.6) 


60.1(7.3) 


56.8(2.7) 


70.8(2.6) 



c 3 [GeV 



QCDSF-UKQCD@ ETM[g(J 
Fit-A Fit-B 
-3.4 -4.7 -3.45 



mo [GeV] 
ci [GeV" 1 ] 
ei(lGeV) [GeV" 3 ] 



0.89(6) 
-0.93(10) 
2.8(1.4) 



0.76(6) 
-1.25(10) 
1.7(0.5) 



0.887(67) 
-1.22(20) 



5. Results at the physical point 

We choose m n , and niQ as physical inputs to determine the up-down and the strange quark 
masses and the lattice cut-off. The SU(2) ChPT analyses are applied to the quark mass dependences 
of m n ,f n and fx taking account of the finite size corrections evaluated at the one-loop level. We 
assume a simple linear quark mass dependences for m 2 K . These points are already explained in 



Sec. |3.2[ We also employ linear forms for the chiral fits of the vector and the baryon masses. It 



is difficult to rely on the heavy meson effective theory [)41[] and the heavy baryon ChPT because of 
their bad convergence behaviors in the chiral expansions [42]. Figure || compares the light hadron 
spectrum at the physical point with the experimental values. The results are encouraging. The 
largest discrepancy is at most 3%, which may be explained by possible 0((£?Aqcd) 2 ) cutoff errors. 
What we are left with is a proper analysis of the p meson and the A baryon as resonance states. 
We calculate the bare quark masses using the axial vector Ward-Takahashi identity (AWI) 



with the nonperturbatively 0(a) -improved axial vector current[]24[]. The renormalization factors 
Za and Zp are perturbatively evaluated up to one-loop level[^, |2q ]. In Table ^ we present the 
results for the quark masses, the lattice cutoff and the pseudoscalar meson decay constants together 
with the experimental values and the recent 2+ 1 flavor lattice QCD results. The physical up-down 
quark mass is just 30% smaller than our lightest one mj^ s (/i = I /a ) = 3.5 MeV at (k" U( j,K" s ) = 
(0.13781,0.13640). We find that the RBC/UKQCD results for the physical quark masses give 
rather larger values than ours. This may be caused by different renormalization procedures: a 
nonperturbative implementation of the RI-MOM scheme for the former and one-loop perturbation 
for the latter. In fact, the ratio of the up-down quark mass to the strange one, which should be 
independent of the renormalization factors, shows a good agreement. It may be also suggestive that 
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Figure 6: Light hadron spectrum extrapolated to the physical point (red circles) in comparison with the 
experimental values (black bars). 



the MILC results obtained with two-loop renormalization factors are found in the middle of ours 
and the RBC/UKQCD ones. We are now evaluating the nonperturbative renormalization factors 
employing the Schrodinger functional scheme. 

In Table [| we also present the results for the pseudoscalar meson decay constants at the 
physical point employing the perturbative renormalization factor at one-loop level. We find that 
our results are consistent with the experimental values within the errors. On the other hand, the 
RBC/UKQCD results are smaller than the experimental values by 5%, albeit rather large systematic 
uncertainties are estimated. Since the MILC results are free from the uncertainties due to the finite 
renormalization, they use f K as a physical input. 



Table 4: Cutoff, renormalized quark masses, pseudoscalar meson decay constants determined with ra,, niK 
and ma as physical inputs. Quark masses are renormalized at 2 GeV. 



physical point 
w/o FSE w/ FSE 



experiment[[43|] RBC/UKQCD [@ MILCQ 



[GeV] 
mjjjP [MeV] 
[MeV] 
m s /m ud 
fn [MeV] 
h [MeV] 
Ik/ fn 



2.176(31) 
2.509(46) 
72.74(78) 

29.0(4) 
132.6(4.5) 
159.2(3.2) 
1.201(22) 



2.176(31) 
2.527(47) 
72.72(78) 

28.8(4) 
134.0(4.2) 
159.4(3.1) 
1.189(20) 



130.7±0.1±0.36 
159.8 ±1.4 ±0.44 
1.223(12) 



1.729(28) 
3.72(16)(33)(18) 
107.3(4.4)(9.7)(4.9) 
28.8(0.4)(1.6) 
124.1(3.6)(6.9) 
149.6(3.6X6.3) 
1.205(18)(62) 



continuum 

3.2(0)(1)(2)(0) 
88(0)(3)(4)(0) 

27.2(1)(3)(0)(0) 
input 

156.5(0.4)(^°) 
1.197(3)(+f 3 ) 



6. Results of the physical point simulation 



A direct simulation at the physical point is attempted by choosing the hopping parameters as 
(fc ud ,fc s ) = (0.137785,0.13660), which is estimated using the four data set with K ud > 0.13754. 
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Figure 7: Chiral extrapolation of m^/m^ 1 together with the result at (fCud,^ 
SU(2) ChPT fit is applied to four data points with m^f^O.Ol. 



(0.137785,0.13660). 



The measured up-down and the strange quark masses are 

m*f (n = 2GeV) = 3.59(20)MeV, mf*^ = 2GeV) = 73.72(44)MeV, (6.1) 

where a~ l is determined from ma at (k u< i, K" s ) = (0.137785,0.13660). The strange quark mass is 
successfully tuned at the physical value given in Table ||. Although the up-down quark mass is 
5a away from the physical value, it is shown that the chiral corrections due to the deviation are 
negligibly small or comparable to the current statistical errors. We discuss the importance of the 
physical point simulation from a view point of controlling the systematic errors. 

We first investigate the pseudoscalar meson sector. In Fig. |7| we plot m\/m^ 1 as a func- 
tion of m 2 n /m^ 1 together with the fit results with the SU(2) ChPT The result at (K ud ,K s ) = 
(0.137785,0.13660) shows an agreement with the extrapolated value at the physical point (red 
star) within rather large statistical error. Figure || illustrates an importance of the physical point 
simulation. For f n we observe that the extrapolated values at the physical point with various fit 
formulae are consistent within the error bars. The story is different for The extrapolated values 
are scattered beyond the error bars. This difference is caused by the typical magnitude of the error 
bar: 4% for f n and 1 — 2% for The increase of the resolution with the diminishing statistical er- 
rors reveals the discrepancies between the various chiral extrapolations. Figure |8| tells us that these 
uncertainties should be removed by a direct simulation at the physical point. In the pseudoscalar 
meson sector it is clear that the up-down quark mass at (k" U( j, fc s ) = (0.137785,0.13660) is so close 
to the physical value that the chiral corrections from the physical point are not relevant within the 
current statistical errors. 

Since we employ mjj as a physical input to determine the lattice cutoff, it is intriguing to check 
the result of niQ at (jc u d, K" s ) = (0.137785,0.13660). In Fig. ^ we compare it with the extrapolated 
value at the physical point. We find a good agreement between them. Moreover, the result at 
(K"ud)K"s) = (0.137785,0.13660) has a smaller error than the extrapolated value. In this case the 
chiral corrections from the the physical point could be comparable to the statistical error. It is 
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Figure 8: Chiral extrapolations of f K (left) and fx (right) together with the results at (rC u d,K: s ) = 
(0. 137785 ,0.1 3660). Red symbols represent the SU(2) ChPT fit applied to four data points with m^f 1 <0.01 . 
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Figure 9: Chiral extrapolation of mji together with the result at (ff ur j, K s ) — (0.137785,0.13660). Linear fit 
is applied to four data points with w^^O.Ol. 



fascinating that we are allowed to make a precise measurement of the physical input directly at the 
physical point. 

For the vector channel we compare the results for m^* and tha at (/c ur j, fc s ) = (0.137785,0.13660) 
with the extrapolated values in Fig. 10 The situations are similar to the niQ case: The extrapo- 
lated values are confirmed by the direct measurements at (/c U( j, ff s ) = (0.137785,0.13660) with the 
smaller statistical errors. 
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